Non-classical correlations in non-Markovian continuous variable systems 
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We consider two identical and non-interacting harmonic oscillators coupled to either two independent bosonic 
baths or to a common bosonic bath. Under the only assumption of weak coupling, we analyze in details the non- 
Markovian short time-scale evolution of intensity correlations, entanglement and quantum discord for initial 
two-mode squeezed-thermal vacuum states. In the independent reservoirs case we observe the detrimental 
effect of the environment for all these quantities and we establish a hierarchy for their robustness against the 
environmental noise. In the common reservoir case, for initial uncorrelated states, we find that only quantum 
discord can be created via interaction with the bath, while entanglement and sub shot noise intensity correlations 
remain absent. 
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I. INTRODUCTION 

Quantum correlations have been the subject of intensive 
studies in the last two decades, mainly due to the general be- 
lief that they are a fundamental resource for quantum infor- 
mation processing tasks. Perhaps, the first rigorous attempt 
to address the classification of quantum coiTelation from an 
information viewpoint has been put forward by Werner |1], 
who introduced an operational definition of quantum entan- 
glement as the property of states that cannot be prepared by 
local operations and classical communication between the two 
parties. One might have thought that such classical informa- 
tion exchange could not bring any quantum character to the 
correlations in the state. In this sense separability has often 
been regarded as a synonymous of classical! ty of correlations. 
However, it has been shown recently |2, 3] that this is not 
the case, and a measure of correlations - quantum discord - 
has been introduced as the mismatch between two quantum 
analogues of classically equivalent expressions of the mutual 
information. For pure entangled states quantum discord co- 
incides with the entropy of entanglement. However, quantum 
discord can be different from zero also for some (mixed) sepa- 
rable state. In other words, classical communication can give 
rise to quantum correlations due to the existence of non or- 
thogonal quantum states. Quantum discord, therefore, cap- 
tures quantum correlations more general than entanglement. 
Separable mixed states having nonzero quantum discord have 
been proven to provide computational speed up in some quan- 
tum algorithms |3, |51] compared to their classical counterpart. 
In addition, the vanishing of quantum discord between two 
systems has been shown to be a requirement for the complete 
positivity of the reduced subsystem dynamics [a]. 

The definition of quantum discord involves an optimization 
problem that, in general, can be tackled only for very simple 
systems. Even in the simplest bipartite system, i.e., a system 
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of two qubits, an analytic expression for the discord for the 
most general two-qubits state does not exist. The optimization 
problem has been indeed only recently solved for the subset 
of states which are unitary locally equivalent to the so-called 
X-states ||3,lal but remains unsolved for more general states. 
For this reason the dynamics of the discord in presence of the 
environment has been up to now studied only for very sim- 
ple finite-dimensional systems. In the case of two qubits, e.g., 
both the Markovian and the non-Markovian time evolution of 
quantum correlations have been investigated [9-12] and it has 
been shown that discord and entanglement behave differently 
under the effect of the environment. In part icular, the phe- 
nomenon of entanglement sudden death fl3], i.e., the com- 
plete loss of entanglement after a finite time, does not occur 
for quantum discord, which instead disappears only asymptot- 
ically |9]. Remarkably, for two-qubits systems in presence of 
nondissipative noise, the discord may remain constant in time 
for very long time intervals, providing the first evidence of a 
quantum property that is completely unaffected by the envi- 
ronmental noise for very long times l\M . The sudden tran- 
sition from classical to quantum decoherence, associated to 
the constant discord phenomenon has been recently observed 
experimentally in a quantum optical set up llisll . 

In this paper we present a detailed analysis of the time evo- 
lution of quantum discord for a more involved bipartite sys- 
tem, namely a system consisting of two non interacting har- 
monic oscillators initially prepared in a thermal twin beam 
(TWB) state. The analytic formula for the quantum discord 
for generic bimodal Gaussian states has been discovered only 
very recently [la, 47]. We use such definition to evaluate how 
the quantum correlations evolve in presence of both indepen- 
dent and common bosonic thermal reservoirs. 

The system we are going to analyze have an immediate ap- 
plication in a quantum optical setting where it may be im- 
plemented by parametric downconversion (PDC), which has 
been addressed as a convenient and feasible setting to visual- 
ize the evolution of quantum correlations 1 18, 19]. In turn, the 
pair of field modes obtained from thermally seeded PDC is as 
a convenient physical system to analyze the quantum-classical 



transition in the continuous variable regime. This scheme 
have been already investigated in ghost- imaging/diffraction 
experiments lUSll . where it has been shown that both entangle- 
ment and intensity correlations may be tuned upon changing 
the intensities of the seeds Ill8l[l91 . In this framework, be- 
sides fundamental quantities like entanglement and quantum 
discord, we will also evaluate a more operational quantity as 
the degree of correlations between the intensities of the two 
beams exiting the noisy channel. The shot-noise limit (SNL) 
in a photodetection process is defined as the lowest level of 
noise that can be achieved by using semiclassical states of 
light [21], that is Glauber coherent states. On the other hand, 
when a noise level below the SNL is observed, we have a 
genuine nonclassical effect. For a two-mode system if one 
measures the photon number of the two beams and evaluates 
the difference photocurrent the SNL is the lower bound to the 
fluctuations that is achievable with classically coherent beams 
and a noise level below the SNL indicates the presence of non- 
classical correlations between the beams. 



II. PHYSICAL MODELS 

In this section we introduce a physical model widely used in 
the description of the non-Markovian dynamics of CV quan- 
tum channels. The main system is made of a pair of identical 
non-interacting harmonic oscillators of frequency ujq and unit 
mass. The free Hamiltonian reads 



Ho = H^, +H°2=\Y. (^i" + ^o^j 



), 



(1) 



At\ 



where Pj = -^ [cij ~ aA and Xj = -7s (oj + aj ) are the mo- 
mentum and position operators, respectively, and aj the field 
operator of the harmonic oscillators (the index j = 1 , 2 labels 
the oscillators). Additionally, we suppose that the harmonic 
oscillators interact with an external environment. In the fol- 
lowing we introduce two different interaction models. 



We consider the case in which system and environment are 
weakly coupled but we do not perform the Markov approxi- 
mation so our results also described the initial short time cor- 
relations between the system and the reservoir. The lifetime 
of such correlations depends on the structure of the environ- 
ment. When the spectral density of the environment changes 
significantly for frequencies close to the system characteris- 
tic frequency, the correlations between system and reservoir 
persist for a longer time and non-Markovian approaches are 
necessary. 

The dynamics of entanglement in such structured reservoirs 
has been studied in both the common 0221 12311 and the inde- 
pendent reservoir scenario [24 l2J. 12511 ■ Here we compare the 
time evolution of the discord with the one of both the entan- 
glement and the intensity correlations. In the case of indepen- 
dent reservoirs, we can establish a hierarchy of nonclassicality 
markers in terms of their robustness against the destructive ac- 
tion of the environment. In the common reservoir scenario we 
find that if the initial state does not possess quantum correla- 
tions, i.e., all three markers of nonclassicality here considered 
have initially zero value, as time passes the interaction with 
the common reservoir can create quantum discord between the 
two system oscillators. Entanglement and nonclassical inten- 
sity correlations, however, cannot be created by the common 
reservoir in the weak coupling limit here studied. Finally, we 
analyze how the quantum discord behaves as a function of the 
initial thermal component of the TWB state. We discover that 
this quantity influences the rate of change of the discord in 
both the independent and the common reservoir cases. 

The paper is structured as follows. In Sec.|ll]we present the 
microscopic physical models for the system and the reservoir. 
In Sec. Hn] we introduce the three markers of nonclassicality 
and, in particular, of nonclassical correlations considered in 
this paper: intensity correlations, entanglement and discord. 
Sec. |IV] investigates the dynamics of the markers in presence 
of common or independent reservoirs. Finally, Sec.lVjcloses 
the paper and draws some concluding remark. 



A. Independent reservoirs 

The first model consists of an external environment made 
of two independent bosonic baths with free Hamiltonian 



Hb = 









"^jkW^kQ^ 



2m 



3 k 



(2) 



The index j = 1,2 labels the bath, and k runs over all the 
bath modes. The Il^fc (Qjk) are the momentum (position) op- 
erators, while Wjk and rrijk are the frequencies and masses 
associated to each bosonic mode. 

Each system oscillator interacts with its own bosonic bath 
(same index j) through a position-position coupling described 
by the following interaction Hamiltonian 



Hi = a^jjkXjQjk, 



(3) 



where -fjk are the coupling constants between the j-th oscil- 
lator and the fc-th mode of its bath and a is a dimensionless 
coupling constant. For the sake of simplicity hitherto we as- 
sume that the baths have the same spectral structure and are 
equally coupled to the oscillators. 

The reduced dynamics of the two oscillators in the case of 
stationary reservoirs is described by the following exact time- 
local master equation lual 



-n{t)[X„[P,,g{t)]] + -r{t)[Xlg(t)] 

lj{t)[X,,{P,,Q{t)}], 



(4) 



where g{t) is the reduced density operator of the oscillators 
and H^ is the free Hamiltonian of the j-th oscillator The time 
dependent coefficients, describing diffusion (A(t), n(i)). 



damping (7(i)) and free frequency renormalization(r(i)) pro- 
cesses, can be expressed as power series in the system- 
reservoir coupling constant a. In the weak coupling limit we 
can stop the expansion to second order in a and obtain ana- 
lytic solutions for the coefficients. We provide their expres- 
sions in the Appendix for high temperature reservoirs char- 
acterized by an Ohmic spectral density with Lorentz-Drude 

cutoff J(w) = ^^^^. 

Using the characteristic function approach the solution of 
(|4|l in the weak coupling limit is given by 112711 

Xt(A) = exp{-A^[W(t) ® W(t)]A} 

xxo(e-^<*^/'[R-'W©R"'(i)]A), (5) 

where Xt(A) is the characteristic function at time t, xo is 
the characteristic function at the initial time i = 0, A = 
{xi,pi,X2,P2)'^ is the two-dimensional phase space variables 
vector, r(i) = 2 J* "f{t')dt', and W(t) and R(t) are 2 x 2 ma- 
trices whose expression is also given in the Appendix. 

The interaction between oscillators and baths is bilinear in 
position and momentum, thus it induces a Gaussian evolution. 
This is of great importance because, as we will see, analytic 
expressions for quantum correlations can be obtained only in 
the case of Gaussian states. 

The characteristic function of a Gaussian state with zero 
mean depends only on the expression of the covariance ma- 
trix a, whose elements are defined as (Tij = ({A.;, Aj})/2 — 

(Ai)(Aj), where A = (Xi, A, -^2,^2) and (•) indicates the 
mean value over the state. 



1 



Xo(A)=exp --AV(0)A 



(6) 



Using Eqs. (|5]l and (|6]l we get the evolution of the covariance 
matrix under the action of the two independent reservoirs 

a{t) =e-r(*)[R(t) ® R(t)]cr(0)[R(t) © R(t)]"^ 

+ 2[w(t)ew(i)]. 

The solution in the case of an initial symmetric covariance 
matrix in its normal form is 



(7) 



a(0) = 



Ao Co 
Co Ao 



a{t)^ 



^md Qmd 



Cind 



, ind 



(8) 



with Aq ~ al, Cq ~ diag(ci,C2), a > and ci, 0-2 real 
numbers, and 1 the 2x2 identity matrix (note that Co and 
CJ."** are symmetric matrices). The analytic expression of the 
matrices A}."'* and C^"^ is given in the Appendix. 



In order to write the reduced dynamics solution g{t), we 
first apply a canonical transformation to the Hamiltonian fol- 
lowing the lines of |[23ll . We define new position X± = 
{Xi ± X2)/V^ and momentum P± = (A ± P2)/V2 op- 
erators for the system. Under this transformation the total 
Hamiltonian becomes 



Ho = 



Pl + P^ 



Hr = 



2 



-{X\ 



Xt 



rukwlQl 



^-^ \2mk 






Hi=aV2X+}_^-fkQk 

k 



(10a) 
(10b) 
(10c) 



In this picture only one oscillator interacts with the bath 
through a position-position coupling, the other evolving 
freely. It follows that the master equation for the reduced state 
g{t) in the new picture becomes 

m ^^[{Hl + Hi), g{t)] - V2A(t)[l+, [X+, g{t)]] 
in 

+ V2Il{t)[X+, [P+, g{t)]] + ^^(*)[^+' ^(*)] 



tV2j{t)[X+,{P+,git)}] 



(11) 



This is of the same form of (01, except for the fact that only 
one effective oscillator is coupled to the environment. The 
dynamics in terms of the characteristic function is then 

Xt(A±) =exp{-A^[V2W(t)®0]A±} 

xxo([e"^R"'W©R"'W]A±), (12) 

with A± = {x+,p^,x-,p-)'^ and being the 2x2 zero 
matrix. Equivalently, the associated covariance matrix d-{t) 
evolves as 



a{t) = e ^/^R(i)®R(t) ct(0) 
+ 2%/2[W(i)®0]. 



V2 



R(t) e R(t) 



(13) 



As in the previous case an initial Gaussian state will maintain 
its character during the time evolution. Indeed, the canonical 
transformation, its inverse and the dynamical evolution are all 
Gaussian operations. 

Given the initial covariance matrix o'(O), applying the trans- 
formations and using Eqs. (fT2l i and ( fTST l we get 



B. Common reservoir 



a(0) = 



Ao Co 
Co Ao 



a{t) = 



A com i^com 

* * 1 (14) 

/^com A com / 5 V^^/ 



In the second example of system-environment interaction 
model, we consider a common bosonic bath, and look at the 
case in which both system oscillators interact with it symmet- 
rically. The Hamiltonian reads 



Hf 



H,= 



E 



\2mk 2 



3,k 



(9) 



with C^"™ symmetric matrix. Details of the solution are still 
given in the Appendix. 



III. NON-CLASSICAL CORRELATIONS 

In the last decades there has been a growing interest in the 
issue of identifying and possibly quantifying the quantumness 



of states of a given physical system. One of the reasons is 
that states possessing quantum features may be useful for cer- 
tain quantum information and computation protocols, or in the 
field of precision measurements, enhancing computation and 
measurements efficiencies. 

In the case of bipartite (or in general multipartite) systems 
the interest is directed not only towards the quantumness of 
the state itself, but also towards the quantumness of correla- 
tions between the different parts. In this paper we provide 
new insight on this issue by comparing different markers of 
quantumness of states and correlations. In particular we are 
interested in studying how non-Markovian dynamical evolu- 
tions affect these quantities in the context of two mode contin- 
uous variable systems. In the following we introduce the well- 
known concepts of intensity correlations and entanglement in 
CV systems as well as the recently introduced quantum dis- 
cord for Gaussian states. 



A. Intensity correlations 

Firstly we consider the intensity correlations marker Icon, 
which is related to the measurement of the two light beams 
intensities, i.e., (fii) and (^2), with fii = {Xf + Pf)/^ being 
the number operator of the i-th mode, and thus feasible with 
current technology. More precisely, the intensity correlations 
marker is defined as lUSl [2811 



{fii + ^2) 



(15) 



For a Gaussian state with zero mean value. Icon depends 
only on the corresponding covariance matrix a and, in the case 
of symmetric Gaussian states, reads 



2fT?3 - al^ - 0-23 - <jI^ 



cm 



(722 — 1 



(20) 



where aij are the covariance matrix entries. 



B. Entanglement 

Entanglement dynamics in dissipative bipartite continuous 
variable domain has been object of interest and numerous 
studies in recent years 11221 - 12511 . Though there exist separabil- 
ity criteria and entanglement measures for a bipartite Gaussian 
state Q (see, e.g., Il30l - [32ll ). in this paper we study the entan- 
glement dynamics by focusing on the logarithmic negativity 
defined as iS 



A/'(e) = max{0, - log(2i>_)}, 



(21) 



with £/_ being the minimum symplectic eigenvalue of the 
partially transpose (PT) covariance matrix of the system, 
namely, a^^ = AaA with A = diag(l, -1, 1, -1). It is 
worth to note that JV{g) > if and only if i^_ < 1/2, that 
is if and only if the state g is entangled: of course, the con- 
dition z/_ < 1/2 is a necessary and sufficient condition for a 
bipartite Gaussian state to be non separable II30I1 . 



and is based on the measurement of the operator I- = fii — 
^2, that is the difference between the intensities of the two 
light modes, whose variance (A/^) may also be written as 

(aP) - (P) - {L)\ (16) 

= (An2) + (An2) - 2(fii)(n2) g^'-\n^,n2). (17) 
where we introduced the second-order correlation function 



g'2'(ni,n2) 



\n1n2) 
("1X^2) 



1. 



In the case of products of coherent states we have Icon — 0, 
which defines the shot-noise limit (SNL) for this particular 
detection process, that is the lowest level of noise that can 
be obtained by using the semiclassical states of light, i.e., the 
coherent states. On the other hand, when 



C. Quantum Discord 

The total amount of correlations in a bipartite quantum sys- 
tem having density operator g is quantified by the quantum 
version of the mutual information 



I{g) = S{g^) + S{g2)~S{g), 



(22) 



where S{-) is the Von Neumann entropy, and gi{2) = 
/I o\ Tr2(i) [g]. Usually the total correlations are divided in a quan- 
tum part, known as quantum discord T>{g), and a classical part 
C{g). The classical correlations are defined as the maximum 
amount of information we can gain on onepart of the system 
by locally measuring the other subsystem ylj, 



C(g) = max 

Hi 



{S{gi)-^P,S{g^^^)], (23) 



< Icon- < 1, 



(19) 



the fluctuations on the intensity correlations are below the 
SNL, indicating genuine non-classical features in the state of 
the system. It is worth stressing that intensity correlations 
below the SNL can be observed also for product states, for 
example in the presence of local squeezing. Hence, this fea- 
ture is not necessary related to the quantumness of correlations 
among different parts of our bipartite system, but rather to the 
quantumness of the overall state itself 1.29.1 . 



where g^L = Tr2(£'l ® 11.;) is the post measurement state 
in which system 1 is left when the result i occurs in a mea- 
surement of system 2 with probability pi = Tri 2(£>ll ® Hi). 
The maximum is taken over the all positive operator valued 
measures {11;} (POVM), X^i n* = ^ performable on one 
subsystem. Classical correlations are thus obtained in cor- 
respondence of the POVM that minimizes the conditional en- 
tropy X^iK'^Cfi'ii^) '■^■' ^^^ allows one to obtain the highest 
amount of information on the state of system 1. The above 



definition is in general non symmetric with respect to the in- 
terchange of the subsystems. In our case however, due to our 
specific choice of the system's initial states (see below) and to 
the symmetry of the coupling with the bath, (|23] | turns out to 
be symmetric during the entire evolution; therefore no specific 
indication of the subsystem measured will be needed. The 
quantum discord is then defined as the difference between the 
total correlations and the classical correlations 



v{e)^i{e)-C{Q). 



(24) 



A peculiar property of quantum discord is that it can be 
nonzero even if the state is separable. This is an indica- 
tion of the fact that entanglement is not the only source of 
quantum correlations. Recently examples of quantum com- 
putational algorithms showing a speedup with respect to the 
classical counterparts, also in the absence of entanglement, 
have been presented IJ, III]. It is believed that the presence 
of quantum correlations other-than-entanglement is responsi- 
ble for this feature. In this sense it is important to study how 
the quantum discord evolves in presence of the external envi- 
ronments, comparing, e.g., its behavior with the behavior of 
entanglement. In the following we answer to this question for 
Gaussian states of CV systems. 

To evaluate the total quantum correlations we use a recently 
developed expression valid only for Gaussian states [16, 17]. 
Given the block form of the covariance matrix a{t) (HI and 
(O, in the symmetric case, the Gaussian quantum discord is 
defined as 

V{g) = f[^ del At) +f(VA^)~f{n+)-f{n^), (25) 

where f{x) — {x+ i)ln(a;+ i) — (x— i)ln(a;— ^), n± are 
the symplectic eigenvalues of the covariance matrix and 



r = At - Ct(At 



CTM 



-^Ct^ 



(26) 



is the covariance matrix of the state of the system A after 
the generalized Gaussian measurement on the system B, de- 
scribed by the covariance matrix 



<:^M 



cosh 2p ^ 1 + tanh 2p cos cj) — tanh 2p sin 
— tanh 2/9 sin </) 1 — tanh 2p cos ( 



(27) 
with p > and < (/i < 2ti. For a generic Gaussian state we 
must perform a minimization procedure in order to find the 
appropriate generalized measurement for the calculus of the 
quantum discord. In the case of Ct = diag(c, — c) the min- 
imum is obtained for a completely heterodyne measurement 
116.1 . For a generic covariance matrix in its normal form the 
exact expression for the discord has been evaluated in lUTJl . 
However in our cases the time evolution of the covariance 
matrix is not in the normal form. Therefore to evaluate the 
discord it was faster to implement a numerical minimization 
procedure. 



IV. RESULTS 

In this section we study the evolution of the previously 
introduced markers of nonclassicality under the influence of 



either independent or common reservoirs, in the weak cou- 
pling regime. We limit our investigation to the case of an 
Ohmic spectrum with Lorentz-Drude cut-off at high tempera- 
tures focusing on the non-Markovian short time-scale. More- 
over we fix hitherto the coupling constant a — 0.1 and 
ksT/hwc — 100, according respectively to the weak coupling 
and high temperature assumptions. 

An important parameter in our discussion is the ratio be- 
tween the cut-off frequency of the baths spectrum ojc and 
the oscillator frequency wq, namely the resonance parame- 
ter X = oJc/wq. In i25ll we have studied the non-Markovian 
entanglement dynamics noting the existence of two dynam- 
ical regimes (x <C 1 and x ^ I) characterized by qualita- 
tively and quantitatively different dynamical behaviors. Gen- 
uine non-Markovian effects occur in the a; <C 1 regime be- 
cause the time-dependent coefficients in the Master equation 
attain negative values in certain time intervals. This feature 
leads to entanglement oscillations, which are not present in 
the X ^ 1 case. As we will see, the same conclusion is valid 
for the intensity correlations and the quantum discord, also 
in the common reservoir scenario. Moreover, in general, the 
dynamics for a; <C 1 is much slower In the following we 
concentrate especially on the x ^ 1 regime, or linear spec- 
trum regime, unless qualitatively different phenomena can be 
reported in the other regime. 

Let us consider as initial states the thermal TWB states de- 
fined as 



g,nir, Ni,N2) = ^2(r)i^(A^i) ^ HN2)Slir) 



(28) 



where i/(iV) = X;„ ^"(l+^)"^"+^V>("-l is a thermal state 

with N average photons and 82(1') — exp{r{a\a\ — 0.102)} 
is the two-mode squeezing operator. In the symmetric case 
we are interested in the two thermal states are characterized 
by the same temperature parameter Ni — N2 — N, while the 
squeezing parameter r can assume any non negative value. 
The covariance matrix a{0) of the initial state is given by 



Ao 



a 
a 



Co 



c 
-c 



(29) 



with a= {N + 1/2) cosh(2r) and c = {N + 1/2) sinh(2r). 
If r = the state is initially uncorrected. If r > the state 
may be entangled or separable, depending on the value of N, 
but it will always possess non-zero quantum discord itlol . 



A. Independent reservoirs 

Initially uncorrected states cannot become correlated at a 
later time when evolving under local operations, as in the case 
of the independent reservoirs model. Thus we focus here on 
initially correlated states (r > 0). In all the various examples 
we examined, we observed that the interaction with the reser- 
voirs has a detrimental effect for all the quantumness markers 
introduced in the previous section. In this sense, not only the 
state becomes more classical but also the quantum correla- 
tions decrease. 



Entanglement however behaves differently from quantum 
discord. Indeed entanglement can disappear after a finite time, 
exhibiting a sudden death and, depending on bath parameters 
and temperature, also exhibit partial revivals |25]. On the con- 
trary, in our system the quantum discord vanishes only for 
t — > oo, a result which is independent from the value of the 
resonance parameter and, at least in the weak coupling limit, it 
is also independent from the spectral distribution and temper- 
ature regime. This is a consequence of the fact that Gaussian 
quantum discord is zero if and only if the Gaussian state is a 
product state and therefore if and only if the determinant of 
the C matrix is zero Ill6lll7ll . This condition is never satisfied 
in weak coupling case for initial two-mode squeezed thermal 
states (see Appendix). 

In FIG. [H we show the behavior of Icon (only the sub-shot 
noise regime), J\f{g) and ^{g) for initial states with r = 2 
and A^ = as a function of cuct. Note that, in the figures we 
have scaled all these quantities so that their initial value co- 
incides. In FIG. [U (top) we choose x — 10, while in FIG. [U 
(bottom) we have x — 0.2. The dynamics is faster when x 
is large and it does not present oscillations typical of a non- 
Markovian evolution. Moreover entanglement is more robust 
to the detrimental effect of the environment than Icorr- This is 
valid also for a; <C 1, a fact that helps us to set up a hierarchy 
for the behavior of our quantities under the influence of the in- 
dependent environment. We will see that such a classification 
cannot be done instead in the more complicated dynamics due 
to a common reservoir 

On the other hand when a; <C 1 non-Markovian oscillations 
are present in all quantum markers. Moreover, in this case, 
intensity correlations and entanglement go to zero at the same 
time, independently from the initial thermal squeezed state. 

We conclude this section with an analysis of the dynamics 
as a function of the thermal parameter N of the initial state, 
concentrating in particular on the dynamics of the quantum 
discord. In FIG. |2] we show the time evolution of 2?(£i) for 
X = 10, r = 2 and different values of A^ = 0, 1, 5, 10. One 
can clearly see that the higher is A^ the slower is the loss rate of 
quantum correlations under the action of local bosonic baths. 
In other words, states with initially higher thermal component 
loose quantum correlations, as measured by the discord, more 
slowly than states with smaller thermal component. 



B. Common reservoir 



N(p) 




FIG. 1: (Colors online) Evolution of intensity correlations marker 
Icorr below the SNL (blue dashed line), logarithmic negativity J^{g) 
(red dashed line) and quantum discord ©(f?) (solid yellow line) in 
the independent reservoir case as a function of ujct. Parameters: a = 
0.1, k.BT/HiJc = 100 and (top) s = 10, r = 2, A = 0; (bottom) 
X = 0.2, r = 0.5, iV = 0. 



D(P) 



0.5 ro 1.5 2.0 ^^4. 

FIG. 2: (Colors online) Evolution of quantum discord ^{q) as a 
function of the scaled time LUct in the independent reservoir scenario. 
Parameters: a = 0.1, ksT/hLj^ = 100, a; = 10 and r = 2. Differ- 
ent lines represent different values of the thermal parameter: A'^ = 
(Solid blue line), iV = 1 (Dashed red line), A = 5 (Dotted-dashed 
yellow line) and A'^ = 10 (Dotted green line). 



In the common reservoir case the CV system dynamics is 
much richer than in the independent reservoir case. Initially 
uncorrelated states, e.g, become in general correlated as time 
passes. For the class of initial Gaussian states considered in 
the paper, however, and in the weak coupling limit, we find 
that, if r = 0, neither entanglement nor intensity correlations 
below the SNL are created (always Icorr < ). On the con- 
trary, quantum discord is created by the action of the common 
environment and grows as time passes, for any the value of 
A^ and of the resonance parameter x, as shown in FIG. [3] for 
X = 10. This result holds also in the Markovian case Il22ll . 
Therefore, in CV systems, the common reservoir always cre- 



ates quantum correlations in the weak coupling regime. 

We note that the initial value of A^ affects the rate of change 
of the quantum discord similarly to the independent reservoir 
scenario discussed in the previous section. While there, for 
initial correlated states, the higher was A^ the slower was the 
decrease of the discord, here, for initially uncorrelated states, 
the higher is A^ the slower is the increase of the discord, as 
one can see from FIG.[3j in the a; = 10 case. This result holds 
for any value of the resonance parameter x. 

We now turn to the initially correlated case, i.e., initial 
states having r > and A^ > 0. In the common reservoir see- 




FIG. 3: (Colors online) Evolution of quantum discord ^{q) as a 
function of the scaled time ojct in the common reservoir scenario. 
Parameters: a = 0.1, ksT/hwc ~ 100, x = 10 and r = (initially 
uncorrelated state). Different lines represent different values of the 
thermal parameter: A'^ = (Solid blue line), A*' = 0.05 (Dashed red 
line), A'^ = 0.1 (Dotted-dashed yellow line). 



nario a comparison between the dynamics of the three mark- 
ers, as the one presented for independent reservoirs, does not 
provide interesting information. The reason is that it is not 
possible to identify a general hierarchy for the most robust 
quantities under the action of the environment. Indeed the 
dynamics is more strongly dependent on the initial state and 
reservoir parameters and the robustness of each of the mark- 
ers changes case by case without exhibiting a general trend. 
Therefore we will present in the following the most interest- 
ing dynamical features of each quantity separately. 

Initial states possessing intensity correlations initially be- 
low the SNL, for small values of N, always lose them com- 
pletely in a finite time. For large x (FIG. E] (top)) and A^ = 
there are no revivals as expected in this dynamical regime. A 
surprising result is, however, that the larger is the initial two- 
mode squeezing the faster the SNL is reached. One would ex- 
pect, indeed, that initial states with higher values of intensity 
correlations initially below the SNL maintain this quantum 
property for longer times than initial states having smaller val- 
ues of intensity correlations initially below the SNL. An oppo- 
site result is reached in the case of small x (FIG. |4] (bottom)), 
where the environment leads to a faster loss of intensity corre- 
lations when r is smaller. In this case however the SNL limit 
is reached at the same time for each value of r and some os- 
cillations and revivals are present due to the non-secular terms 
of the solution. 

The behavior of entanglement in these systems has been 
studied in detail in previous papers ll22ll23ll . We can however 
summarize the most important features. In the weak coupling 
regime, if there is no initial entanglement, it is not possible 
to create it. This has been already pointed out previously for 
uncorrelated initial conditions. The same conclusion holds, 
however, also for initially correlated states with high enough 
values of N in the initial state. In general, for initially en- 
tangled states, in the short-time scale we can observe sudden 
death and revivals depending on the values of A^, r ,x. If the 
initial value of the entanglement is small, it is rather difficult 
to observe revivals. Some revivals can be seen in the x <C 1 
case, usually due to non-secular terms |25]. When the sys- 
tem is initially strongly correlated entanglement sudden death 




FIG. 4: (Colors online) Dynamics of the intensity correlations 
marker Icorr for a common reservoir as a function of ujct. Param- 
eters: a = 0.1, kBT/huc = 100, N = 0, (top) x = 10; (bottom) 
X = 0.3. r = 0.5 (Blue solid line), r = 1 (Purple dashed line), 
r — 1.5 (Yellow dot-dashed line), r — 5 (Green dotted line). 



and revivals can be observed as well as situations in which in 
the non-Markovian time scale the entanglement never goes to 
zero. 

Finally, we consider the behavior of the quantum discord. 
As in the uncorrelated case, the value of A^ in the initial state 
influence the rate of discord changes. High values of A^ make 
the correlations more robust to the influence of the reservoir. 
The value of the resonance parameter x does not influence the 
qualitative behavior of quantum discord. Small x lead to a 
slower dynamics and presence of oscillations in the solution. 
So we consider as an example the case A^ = and x = 10 
(See. FIG. |5]). The different curves correspond to different 
values of r. If r is small (also r = 0) the discord has an 
initial small value and starts to increase during the dynamics 
for short times. On the contrary if r is large the discord has 
an initial high value and it decreases in time. However, in- 
dependently from the value of r, after an initial transient time 
interval, all the discord curves tend to overlap. The interaction 
with the common bath seems to destroy the original informa- 
tion on the initial correlations during such transient time inter- 
val, driving the CV system towards states which have similar 
value of the quantum discord. 



V. CONCLUSIONS 

In this paper we have studied the time evolution of three 
different indicators of quantum correlations for a bipartite CV 
system initially prepared in a thermal TWB state and interact- 



8 



D(P) 



1.0 




FIG. 5: (Colors online) Evolution of quantum discord ^{q) as a 
function of the scaled time ujct in the common reservoir scenario. 
Parameters: a = 0.1, ksT/huJo = 100, a; = 10, TV = 0. Different 
lines represent different values of the squeezing parameter: r — 
(Solid blue hne), r = 0.2 (Dashed red line), r = 0.4 (Dotted-dashed 
yellow line) and r — 1 (Dotted green line). 



ing with either independent reservoirs or a common reservoir 
at high temperature. Using a definition for the quantum dis- 
cord in CV systems recently introduced in Refs jldllTll . we 
have calculated analytically the dynamics of such quantity un- 
der the only assumption of weak coupling between the system 
and the environment. Moreover, we have compared the dy- 
namics of the discord with the dynamics of entanglement and 
intensity correlations. 

We have demonstrated that, in the independent reservoir 
scenario, initially correlated thermal TWB states, more specif- 
ically states with nonzero initial discord, loose their quantum 
correlations slower and slower for increasing values of their 
initial thermal component A^ . Similarly, in the common reser- 
voir scenario, for initially uncorrected states having zero dis- 
cord, the higher is N, the slower is the reservoir-mediated rate 
of increase of quantum discord. 

A comparison between the dynamics of discord, entangle- 
ment and intensity correlations shows that, when the two sys- 
tem oscillators interact with independent reservoirs, initial 
nonclassical intensity coiTelations disappear faster than the 
initial entanglement, which in turn disappear faster than the 
initial discord, for x ^ 1. In the opposite regime a; <C 1, the 
dynamics of the intensity correlations and of the entanglement 
is very similar and both quantities disappear at the same finite 
time. 

In the common reservoir scenario, we have studied whether 
initial states possessing zero correlations can become corre- 
lated via the action of the reservoir in the weak coupling 
case here considered. Our results show that only discord 
can be created, while entanglement and initial correlations re- 
main zero if they are initially zero. In fact, quantum states 
with nonzero discord are much more common than entangled 
states, as demonstrated, e.g., in II34I1 . In this sense we expect 
quantum discord to be easier to generate than entanglement. 

Our results are the first attempt to characterize the time evo- 
lution of discord and intensity correlations, comparing them to 
the time evolution of the entanglement, in CV non-Markovian 
systems. Hence they provide a first step in the description 
of the behavior of quantum correlations and their robustness 
under the effect of dissipative environments. In the future we 



plan to investigate whether effects such as the constant discord 
and the sudden transition form classical to quantum decoher- 
ence 13, recently discovered in bipartite qubits systems, may 
occur also in CV systems. 
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Appendix: The master equation coefficients 

The time-dependent coefficients of the master equations ^ 
and ( ITTT i. in the case of thermal reservoirs and at the second 
order in the coupling a are given by 

A{t) ^o?j ds dujJ{Lj)[2N{Lj) + 1] cos(ws) cos(a;os), 
Jo Jo 

(A. la) 

n{t) = a^ / ds dojJ{uj)[2N{uj) + 1] cos(u;s) sin(wos), 
Jo Jo 

(A. lb) 

j{t) — a ds dujJ{uj)sm{ujs)sm{ujos), (A.lc) 

Jo Jo 

l*t /'+00 

,2 



r(t) — a I ds dujJ{uj)sm{ujs)cos{uJos), (A. Id) 

Jo Jo 

with N{uj) = [exp{haj/kBT) — 1]~^ being the mean number 
of photons with frequency oj, while J{uj) defines the spectral 
distribution of the environments. For the Ohmic distribution 
with Lorentz-Drude cut-off in the high-T limit (2N{uj) + 1 ~ 
ksT/huf) we have 

'^\^) ~ nn — T\~ — t'^ ~ ^ [a;cos(T/a;) — sm(r/a;j}, 



(A.2a) 



n(f) 



2(l + a;2) huJc 

— -— ^- {1-e cos(r/j;)+a;sm(r/2;) }, 

2(1 + x'^) huJc 

(A.2b) 



2 2 

lit) = 2(1 + x'^) ^^ ~ e'^[coa{T/x) + 2;sin(r/a;)]}. 

(A.2c) 

We do not provide the analytic expression of r(t) because its 
contribution to the solution in the weak coupling limit is neg- 
ligible. 



Appendix: The master equation solution 

The time evolution of the characteristic functions in our two 
decoherence models (|5]i and ( fT2] ) contain the 2x2 matrices 



R(t) and W(i). Under the weak coupling assumption their 
expressions is 



R(i) 



cosojo^ sinwo^ 
— sin Wot cosuiot 



,-r(t) 



R(t) 



=-r(^ 



'M{s)ds 



(A.l) 



R'^(t), (A.2) 



W{t) = e 
with 

M{t) = R^(i) ( _^g/2 ""?^^' ) ^(^)- (^-3) 
If we expHcit the calculations the following functions appear 

r{t) = 2 ( -/{s)ds, (A.4a) 

Jo 

Ar{t) = e-r(*) /" e^('')A(s)ds, (A.4b) 

Jo 

AcoW==e-^(*' / e^(")A(s)cos[2wo(i-s)]ds, (A.4c) 



A,,(i) 



- p-r(t) 



s^(") A(s) sm[2iOa{t - s)]ds, (A.4d) 



nco(i) = e-^(*) /" e^(^)n(s) cos[2a;o(i - s)]ds, (A.4e) 
n,.i(t) = e^^(*) / e^('''n(s) sin[2wo(i - s)]ds. (A.4f) 

The last five coefficients can be evaluated numerically. How- 
ever if we are interested in the short non-Markovian time 
scale, we can use the approximation exp(±r(t)) ~ 1. Under 
this assumption we can evaluate all the coefficients exactly in 
the case of the Lorentz-Drude spectral function at high-T. 

The last four coefficients are called non-secular. In some 
dynamical regimes their contribution is not essential and can 
be neglected (secular approximation). In this paper however 
we do not perform this approximation. 



2. Common reservoir solution 



To obtain the solution for the common reservoir model we 
first have to transform the original covariance matrix o'(O) into 
the one in the new picture a'(O). Then using Eqs. (fT2l) and ( fT3l ) 
we evolve the matrix into (T{t) and finally we apply the in- 
verse picture transformation to get the solution to the problem 
(T(t). The transformation expressions can be easily obtained 
by comparing the following definitions of the covariance ma- 
trices and implementing the canonical relations 



a., = ({A„A,})/2-(A,)(A,) 



J/' 



a., = ({A±,Af})/2-(A±)(Af), 



(A.6a) 
(A.6b) 



with A = {Xi,Pi,X2,P2) and A± = (X+,P+,X_,P_ 
The solution then reads 



ait) 




(A.7) 



where 



X = ff+(r) a - g-{T) c cos(2a;) + 

y = g+{T) a + g-{T) c cos(2a;) + 

2 = g_ (I j c sin(2xj -= 

V2 



Ar + (Acq - n,,) 

(A.8a) 

Ar - (Aco - Us^) 



V2 



(A.8b) 
(A.8c) 



1. Independent reservoir solution 

As we already pointed out, for an initial Gaussian state with 
zero mean, the solution of the master equation is obtained giv- 
ing the time evolution of the covariance matrix. In the case of 
independent reservoirs we need to apply definitions (|5]i and 
(|6]l to Eq. O. If we do this we find that the covariance matrix 
(O at time i can be written as 



, ind 



-lind 



Aoe" 



A^- 



Arl 



-A,, + n. 



-A, 
-(Ac 



(A.5a) 



sinh(2r) e 



r / cos 2ajo^ sin 2^0^ 



sm2ujot 



cos 2cjoi 

(A.5b) 



and 



fi = -.g-(r) a + .g+(r) c cos(2x) + 



i^ = -ff-(r) a - 5+(r) c cos(2x) + 



Ar- 



(A.9a) 

(Aco - Us^) 



V2 



e = -g+(r)csin(2x) 



A. 



n. 



V2 



(A.9b) 
(A.9c) 



with .g±(r) = 5(1 ± e-^) 
c={N + 1/2) sinh(2r). 



{N + 1/2) cosh(2r) and 
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